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ABSTRACT 



A method is provided for determining Young's modulus and, 
if desired, Poisson's ratio of a coating on a substrate wherein 

load-displacement indentation data in the elastic region 
(either elastic loading or unloading) generated using an 
indenter is analyzed to interpret the elastic repsonse of a 
coated material. The data analysis pursuant to the invention 
is based on an extended Hertzian analysis developed for 
thin-film coatings. 

17 Claims, 11 Drawing Sheets 
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Preparation for measurement, such as 

■ Sample preparation 

■ Instrument calibration 

■ Coating thickness measurement 





Record indenter's 
radius, if available 




■ Profilometer 

■ Scanning electron microscopy 
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Calibration of the indenter 

■ Execute indentation on a 
calibration sample, e.g. a fused 
silica sample 

■ Collect calibration curves of 
load and displacement 

■ Using Eq. (1), determine radius 
with known properties and/or 
vice versa 


Record indenter's 
properties, if available 


< 
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Execute indentation on users' 
samples and collect data points or 
curves of load and displacement 



Minimization 

■ Based on values of each data point of load and 
displacement, apply Eqs. (1) and (2) to find Es* 

■ With a expressed in Eq. (11), minimize Eq. ( 1 2) 
for unknown material properties, which could be 
as many as four: Ec, Vc, Eb and 

■ Average values obtained for these unknowns 
and evaluate Ef„r 

■ Option for Emr < 1 : Use Emr value and Fig. 3 to 
evaluate a. Then repeat the minimization 
process in order to obtain better results 



Output Young's modulus 
and Poisson's ratio 



Fig. 13 
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DETERMINATION OF YOUNG'S MODULUS 
AND POISSON'S RATIO OF COATINGS 
FROM INDENTATION DATA 

This application claims benefits and priority of U.S. 
provisional application Ser. No. 60/511,487 filed Oct. 14, 
2004, the entire disclosure of which is incorporated herein 
by reference. 

CONTRACTUAL ORIGIN OF THE INVENTION 

Funding of the invention was provided in part by the 
OflSce of Naval Research under Grant No. NOOOl 4-01-1 - 
0392/POOOOl. The Government may have certain rights in 
the invention. 

FIELD OF THE INVENTION 

This invention relates to a method of determining Young's 
modulus for coatings with uniform coating tliickness and, 
more particular, to a method of determining Young's modu- 
lus and optionally Pois son's ratio for thin- film coatings 
using load-displacement indentation data in the elastic range 
in an efficient, accurate and non-destructive manner. 

BACKGROUND OF THE INVENTION 

Young's modulus, or elastic modulus (E), and Poisson's 
ratio (v) are two of the key mechanical properties of mate- 
rials. Modem thin-film technologies have provided means 
for the development of coatings for various applications in 
mechanical, electrical, and biomedical engineering to pro- 
tect the surfaces of components that bear load and transmit 
motion for performance and life enhancement. Such coating 
development demands an efficient and reliable, yet conve- 
nient, technology for the determination of the Young's 
modulus and Poisson's ratio of the coatings. 

Currently, five tecluiiques are coromonly used for the 
measurements of the Young's modulus of thin-film coatings: 
(1) indentation, (2) beam bending, (3) vibration, (4) ultra- 
sonic surface wave, and (5) scamiing acoustic microscope 
(References: Bhushan 1999, Fischer-Cripps 2002, Schneider 
et al. 1992, and Fang 1999). Among these methods, inden- 
tation is the most popular technique, with which Young's 
modulus is often determined through the initial part of the 
unloading curve of an indentation data with sharp-pointed 
indenters (References: Loubet et al. 1984, Doemer and Nix 
1986, Pharr et al. 1992, and Oliver and Pharr 1992). How- 
ever, determining Young's modulus of a coating through tliis 
technique, at the macro, micro, or nano scale, has the 
following drawbacks. (1) The indentation data are indis- 
pensably contaminated by the response of the substrate 
materials under the coating to be measured. In order to 
minimize the substrate influence on properties measurement, 
the maximum indentation displacement (i.e. mutual 
approach or depth) is empirically limited to one-tenth of the 
coating thickness. This empirical requirement is not suitable 
for modulus measurement (References: Chudoba et al. 
2002). Empirical and theoretical expressions have been 
proposed to account for the substrate influence on the 
measurement, e.g. References: Doemer and Nix 1986, Gao 
et al. 1992, and Mencik et al. 1997. (2) The indenters have 
to be made sharp in order to have meaningful coating 
response under such shallow indenter penetration. The tip 
shape often deviates from the ideal one (rounded) and the 
indenters are easily damaged. (3) The indentation process is 
destructive. Dents are left in the sample after examination. 
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In the neigliborliood of dents, materials undergo the work- 
hardening process due to plasticity. For ceramic coatings, 
coating cracking usually accompanies the indentation pro- 
cess and affects the result of measurement. (4) For polymers, 

5 creeping may occur in the loading-unloading process and 
affect the result of measurement. (5) Existing models cannot 
accurately determine the contact area. For example, material 
pileup is significant in the indentation of thin-film and 
affected by the residual stress in the sample. The material 

10 pileup and plastic defomiation under the indenter also 
change the coating thickness on the measurement site. 
(References: Lim et al. 1999). (6) Measured modulus is a 
modified Young's modulus, which equals Young's modulus 
divided by the difference between unity and square of 

15 Poisson's ratio. 

Indentation displacement and force could be measured 
very accurately. For example, a commercial nanoindentation 
system could achieve the displacement resolution of 0. 1 nm 
and the force resolution of 1 |liN. The data analysis of 

20 recorded load-displacement indentation data is critical to the 
testing goal, such as determination of Young's modulus. 
Traditional indentation is in elastic -plastic range and 
Young's modulus is determined by elastic punch models 
(References: Loubet et al. 1984, Doemer and Nix 1986, 

25 Pharr et al. 1992, Gao et al. 1992, and Oliver and Pharr 
1992) with the projected contact area and the slope of the 
load-displacement curve at the beginning of the unloading. 
Poisson's ratio for the sample is not determined but esti- 
mated. References (e.g., Wallace and Ilavsky 1998; Herbert 

30 et al. 2001) showed that Hertzian theory could be applied to 
determine the modified Young's modulus of a homogeneous 
sample by a spherical indenter in the elastic range. For 
thin-film coatings. References: Chudoba, Schwarzer and 
Richter (2000) showed the possibility to determine the 

35 modified Young's modulus by a spherical indenter working 
in the elastic range. An analytical solution for the indentation 
displacement due to an estimated (inaccurate) Hertzian 
pressure distribution is utilized to calculate the modified 
Young's modulus of a thin-film through a curve-fitting 

40 process of the load-displacement curve obtained from a 
depth-sensing measurement. However, the assumed Hert- 
zian pressure distribution is not the true contact pressure 
between the spherical indenter and the coated sample, and 
because of the lack of the relationship between displacement 

45 and load a curve-fitting process has to be used. In all known 
indentation data analysis, Poisson's ratio is always assumed 
a priori. 

SUMMARY OF THE INVENTION 

50 

The present invention provides a method for determining 
the Young's modulus and, if desired, Poisson's ratio of a 
coating on a substrate wherein load-displacement indenta- 
tion data generated in the elastic region (either elastic 

55 loading or unloading) is analyzed to interpret the elastic 
response of the coated substrate. The data analysis is based 
on an extended Hertzian analysis developed for thin-film 
coatings pursuant to the invention. 

In one embodiment of the present invention, a method of 

60 the invention involves the steps of generating load versus 
displacement indentation data using an indenter to indent a 
coated sample having a coating on a substrate in an elastic 
range of the coating; determining a modified Young's modu- 
lus of the coated sample (E*^) using the indentation data; 

65 and determining the Young's modulus of the coating (EJ 
using a relationship that expresses the modified Young's 
modulus of the coated sample (E*^) as a ftmction of the 



10 



us 7,165,463 B2 
3 4 

Young's modulus and Poisson's ratio of the coating, the DESCRIPTION OF THE DRAWINGS 

Young's moduhis and Poisson's ratio of the substrate, and a 

parameter related to thickness of the coating. In a particular piG. la shows schematically the indentation testing of a 

embodiment, the relationship is expressed as: coated sample with a flat surface for measuring of the 

^ Young's modulus and Poisson's ratio of the coated sample. 
The coated sample comprises a coating on a base or sub- 
^ _ ^ 1 - + ^ + )g + A/gg^ strate as shown. In FIG. la, the coating with thickness t^, is 

1 + 4aKe - Xk^ deposited on the fixed base and the total thickness of sample 

is t^. The coating has mechanical properties; Young's modu- 
lus E^ and Poisson's ratio v^, while the base or substrate has 
where Young's modulus E^, and Poisson's ratio V^. The indenter 

has a radius of R and its mechanical properties are E^ and v,. 
Before load is applied, the indenter and the sample touch 
Ec{l + vt) each other but there is no indentation displacement (the 

^ ~ Et{i + Vc) dotted lines). After load P is applied, both the indenter and 

the sample are deformed (the deformation of the indenter is 
not shown in all figures). The contact radius is denoted by a 
is the shear modulus ratio* indentation displacement is denoted by h. The 

indentation displacement is the traveling distance of a mate- 
rial point in the indenter far way from the contact interface, 
r-\ 4(1 _ e.g. the center of the sphere. The load and the displacement 

r-\-{?>-Avc)^^ A = 1 - ^ _ 4^^^ ' ai-e recorded at specific steps when the load is increased. 

FIG. lb is a schematic load (P) versus displacement (h) 
curve generated by the above described indentation test. 

25 

and e=exp(-2a) where parameter a is related to non- FIG. 2 shows numerical results of normalized displace- 

dimensional coating thickness. i^ent (h/ho) versus normalized coating thickness (TJ with 

-ru • ^' ^ u ^- A ^ A ^ • different modulus ratios (E^^). 

Ine mvention lurther can be practiced to determme the ^ ^^'^ 

Young's modulus and Poisson's ratio of the substrate using ^ ^^^^^ parameter a for all conditions investigated 

this relationship numerical analyses (FIG. 2). Parameter a is deter- 

. . * ^. . . i . i mined by using displacement data in FIG. 2 and Eq. (Equa- 

In practicing embodiments of the invention, the mdenter ^-^j^^ ^^q^^ 

can include an uncoated surface that indents the sample or t-^^^ a ^ ^^ ^ ^- i • i • ^- x- 

, ^ FIG. 4 shows the representative a, wliich is a lunction oi 

a coated suriace that indents the sample. . . +• i .c n ^- » 

^ dimensionless coatmg thickness, for all ratios of Young s 

The present invention is useful to accurately, conveniently 35 modulus shown, 

and non-destructively determine Young's modulus and Pois- pj^. 5 shows the relative error of predicted displacement 

son's ratio of a coating (including thin-film coating) con- ^sing the extended Hertizian theory in percentage and the a 

currently through indentation testing. The method can be values of FIG. 4 where solid lines are E,„, greater than 1 and 

used as a stand-alone analysis method or can be incorporated dashed lines are E^^ less than 1 . In order to verify that data 

into and practiced in conjunction with existing nanoinden- 40 piG. 4 works for other modulus ratios E^^, Eq. (10) is 

tation systems. used to predict the displacement. The relative error is 

The present invention is advantageous as follows: defined as the difference between predicted and numerical 

(1) During the indentation measurement, curve fitting is ^^^^^^^ ^^^^^^^ numerical resuhs. 

eliminated V\G. 6 shows the configuration of a coated indenter in the 

/^x ^1 • J ^ ^- ^ • J ^ J • ^1 indentation testing. The base of the indenter has mechanical 

(2) The indentation measurement is conducted in the ^ a a ..u ^ ^^ x. a a 

^ \ . IIP - 1 properties En and Vn, and the coating perlectly bonded on it 

elastic range ol loading and therelore is non-destruc- 1 V a -t^. r^i. • j \ 1 ^1 • 1 

^ . -i i 1 r A • has E, and v.-. The coatmg of the mdenter has a thickness of 

tive. No crack, material pileup, or dent ol the coatmg . r^+t,^^ ^oJo«,..+™ c^^^ uir- 1 

' r ° t. Other parameters are the same as in FIG. 1. 

occurs. 

.^.^^ / ,j , , J, FIGS. 7i3f and 7^ show the conventional measurements of 

(3) Noplasticity IS involved, so that the Young s modulus 50 iaija-i- a ^ 

^ \ ^ / . . . i 1 1 Ml 1 wear scars lor coated samples Al and A2 in order to 

haMentn determine coating thickness. FIG. 7c shows schematically a 

ar ening e ec . ^^^^ ^-^j^ radius of 10 ram that is used to remove each 

(4) A preferred method embodiment practiced with coating. After measuring the diameters of the corresponding 
spherical indenters does not require the estimation of circles, one can determine the coating thickness; e.g. sample 
contact area, which usually involves large error. a1 and A2 have a respective coating thickness of 2.08 and 

(5) The indenter is a spherical ball or other body with a 3.15 microns. 

spherical tip, easy to make and can be repeatedly used. piG. 8 shows indentation data of load vs displacement for 

Furthermore, any hard materials with or without a each of coated samples Al and A2 from nine measurements 

coating can be used to make the indenter ball. There- ^sing a sapphire indenter ball of 122 microns radius and E. 

fore, the indenter cost is reduced. of 441 GPa and v, of 0.234 as shown. The data is the average 

(6) The method works for testing any coating, which can of the nine tests in the loading portion with a displacement 
be hard or soft, and comprise metallic, ceramic, or less than 200 nm. Error bars show the standard deviation of 
organic materials, and any substrate material. the nine tests. 

These and other advantages of the present invention will 65 FIG. 9 shows the relation between modified Young's 
become more readily apparent from the following detailed modulus of coated samples Al and A2 and normalized 
description taken with the following drawings. coating thickness (TJ obtained from average indentation 
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data shown in FIG. 8. Results for greater than 1 , which are assumed to be nominally flat. If the sample is free of plastic 

corresponding to measurements with small displacement, deformation, R'^' simply equals R, the radius of the indenter. 

are not shown. Otherwise the radius, R^, of the impression due to plastic 

FIGS. 10a and 10b show predicted Young's modulus deformation will contribute to R*, 

values and Poisson ratio values, respectively, for sample Al 5 
for each discrete indentation measurement. 

FIGS. 11a and 11^ show predicted Young's modulus J_ = i _ J_ 

values and Poisson ratio values, respectively, for sample A2 R* R Rr 
for each discrete indentation measurement 

FIG. 12a shows a generic indentation test using a non- 10 

spherical tip. FIG. 12h shows a corresponding loading- K could be measured from imaging techniques. The 

unloading curve, where the elastic loading may be only a mechanical properties and the radius of the indenter are 

small portion and plastic deformation is unavoidable before often given beforehand or could be calibrated as described 

the elastic unloading. hereafter. With the load-displacement curve, one can deter- 

FIG. 13 is a flow chart illustrating a method aspect of the 15 mine E, , from Eq. (1) and E*, from Eq. (2). By assuming a 

present invention. proper Poisson' s ratio of the sample, the Young's modulus 

of a homogeneous sample can be determined. 

DETAILED DESCRIPTION OF THE However, for coated bodies, the classic Hertzian theory is 

INVENTION applicable to determine mechanical properties of the 

20 coating. 

The present invention provides a method for determining ^he present invention involves extending or modifying 

the young's modulus and, if desfred, Poisson's mtio of a ^he classic Hertzian analysis to enable determination of 

coating on a substrate (piDviding a coated sample) as well as mechanical properties (i.e. Young's modulus and Poisson's 

of the substrate wherein load-displacement indentation data ^^^o) ^^^^ed substrates or bodies (coated samples) lUus- 

in the elastic region (either elastic loading or unloading) ^5 trated in FIG. 1. For example, the coated sample is iUus- 

generated using an indenter is analyzed to interpret the Crated including two different parts: a substrate (base) and a 

elastic response of the coated sample. The data analysis thin-film coating, which are bonded perfectly to each other, 

pursuant to the invention is based on an extended Hertzian ^heir mechanical properties (i.e. Young's modulus E and 

theory developed pursuant to an embodiment of the inven- Poisson's ratio v) are labeled by subscripts b and c for the 

tion for thin film coatings ^^^^ (substrate) or the coating, respectively. The contact 

As background, the classic Hertzian theory presents exact ^ofs^quence of the indenter and the sample should be 

. . 1 n . . j-1 + j between two limitins; cases: the contact ol the indenter and 

contact solutions lor contact pressure, displacement, and , . • i / ^ ^ ^^ ^ • ^ • • i • i 

^ ^ 1 ^ ^ A J ' ^ • n • ^- 1 1 the sample with (a) zero and (b) innnity coatine thickness, 

contact area between two elastic, isotropic, Irictionless sol- ^ ^. . . \- ^ ^ ^ • 

I- ^ ^ 11 1 11 ij- Both limitine cases can be predicted by the Hertzian analy- 

ids in contact and has been widely applied in engineenng .„ ^ .^^ ^-i ^ 

- ^ A- A. A-f^ A \r ' 35 SIS. The contact radius m the first limitmg case (with zero 

practice, includins; measurements ol the modified Youne s . ^.i ii i-i i ii 

11.1 1,5+1. coating thickness) is denoted by ao, which can be calculated 

modulus through depth-sensing tecmiiques, e.g. Relerences: , ^ ^ - i • 

Wallace and Ilavsky (1998). The formula used for determi- ^® Hertzian analysis, 

nation of the modified Young's modulus is the relationship 

between load and displacement derived by Hertz (1882), 



40 (5) 



3P (1) 



HR*h^ 

where E^. ^ is the reduced Young's modulus, consisting of the 
modified modulus of the indenter and the base material, 

where h is the displacement and P is the load applied on the 
indenter with a radius R. E^. ^ is the reduced Young's modulus 
of the indenter and the sample, J_ = J_ + _L 



Eis Ef E; Quantity ag is used to normalize coating thickness, t^. The 

dimensionless coating thickness is expressed as 

55 r,=r/«o (7) 

Subscripts i and s indicate the indenter and the sample, i ^ . ^ i i ^ 

1 .1 • + J + A'f^ A For the contact between the indenter and coated body 

respectively. E with a superscript * denotes the modmed . . i n • . r i i i . .i /• n 

Young's modulus of homogeneous materials, • T """^ l^^t ST- 

Eqs. (1) and (2) onginally presented by the classic Hertzian 

anaylsis, but redefine E*^ in Eq. (2). Thus 



60 



E (3) 

1-V2 



3P 

Eq. (1): = 



R* in Eq. (1 ) is the modified radius. With known R* and E*^-, 65 
a data point on the load-displacement curve gives a E*^ 
value from Eqs. (1) and (2). The sample is generally 



1 1 1 

Eq. (2): — = — + — 
^ Eis EJ E; 
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-continued 

l-(A + /f + W)0 + Aj^ (8) 



where 



is the shear modulus ratio; 



= — :; — : and A = 1 - — — - — - 

r + (3-4vJ l + r(3-4vj 



are similar to the Dundurs parameters (see References: 
Johnson, K. L. 1985); and e=exp(-2a). Here, Equation (8) 
is derived by comparing the normal deformation expression 
for a layered medium in the frequency domain (see Nogi and 
Kato, T., 1997 "Influence of a Hard Surface Layer on the 
Limit of Elastic Contact-Part 1: Analysis Using a Real 
Surface Model", ASME Journal of Tribology, 119, pp. 
493-500) with that for a homogenous medium (the same 
expression, but with the film tliickness being zero). Param- 
eter a is related to the non-dimensional coating thickness, T^. 

Note that E*^ is now a function of coating thickness and 
material properties of the coating and the substrate. The 
essence of Eq. (8) is that the coated sample is equivalent to 
a homogenous virtual sample with a modified modulus 
depending on the dimensionless coating thickness. When a 
is infinity (no substrate at all and 6 is zero), it is evident that 
E*^ in Eq. (8) equals E*^. hi this case, the sample is solely 
composed of the coating material. On the other hand, it can 
be verified that when a is zero, the fraction in Eq. (8) can be 
simplified as. 



i-iX + K) + XK _ e; (9) 
1 -Xk ~ E; 
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to be rigid. Note that in the literature, rigid indenter assump- 
tion is widely used for derivation purpose. A load (P) of 400 
N is applied on the indenter. These parameters detennine 
ao, »0.254 um. The dimensionless thickness of the coating 

5 (TJ varies from 0.01 to 3. Poisson's ratio (vj of the coating 
is set as 0.3 and E^^ (E^=E*JE*fj) varies from 0.25 to 4. 

The corresponding displacement h for the combination of 
conditions selected in the preceding paragraph is evaluated 
using a numerical code or analysis of the type developed by 

10 Liu and Wang in "Studying Contact Stress Fields Caused by 
Surface Tractions with a Discrete Convolution and Fast 
Fourier Transfomi Algorithm", Journal of Tribology, Vol. 
124, pp. 36-45 2002, and in "A Three-Dimensional Ther- 
momechanical Model of Contact Between Non-Conforming 

15 Rough Surfaces", Journal of Tribology, Vol. 123, pp. 17-26 
2001 , the teachings of both of which are incorporated herein 
by reference. The numerical analysis involves an iteration 
process to modify pressure distribution (initial value is 
simply external load divided by simulation area) between 

20 two bodies in contact so that contact conditions could be 
fulfilled. These conditions are: a) the external load should be 
balanced by interfacial force due to the resolved pressure, b) 
when two surfaces of the bodies in contact are separate, 
pressure should be zero; otherwise c) pressure should be 

25 greater or equal to zero. The input parameters are external 
load, geometry and material properties of each body, and 
coating thickness. The numerical analysis determines con- 
tact pressure, displacement h, gap between two surfaces, 
among other outputs. 

FIG. 2 shows the displacement values obtained for the 
above-described combination of conditions, the displacment 
h being normalized by ho, which is the displacement under 
the same conditions except zero T^. If T^ is large enough, it 
is expected that each curve in FIG. 2 should reach a certain 
value. It is apparent that the substrate effect is significant 
even when the coating thickness is three times of Oq. Similar 
to Eq. (1), one can write 

40 3p 

Eib = . 

4y 



45 Therefore, the normalized displacement can be written as 



which equals 1/E^^, where E^^=E*^E*^ is defined as the 
coating-to-base modulus ratio. Thus, in the case of zero a 
(no coating at all), one can find E*^ in Eq. (8) equals E'^^ 
with the help of Eq. (9). Indeed, these two special cases are 
the two limiting cases mentioned above. Equation (8) 
together with Eqs. (1) and (2) pursuant to the invention 
extend the Hertzian analysis to solve indentation problems 
with coated samples. 

The parameter a appearing in 9 and Eq. (8) is determined 
beforehand with numerical simulations described hereafter. 
The extended Hertzian analysis consists of Eqs. (1), (2) and 
(8) and suitable a. In the following, a procedure with 
numerical simulations is presented to determine a before- 
hand. Since numerical simulations need specific conditions, 
an arbitrary set of parameters is chosen in the beginning. 
However, it is shown in the end of the procedure that the 
results of a work well for a wide range of conditions. The 
flat base (substrate) is assumed to be made of steel and has 
Poisson's ratio (v^) of 0.3 and Young's modulus (E^) of 200 
GPa. The indenter has radius (R) of 12 mm and is assumed 



ho \Eis> [E^iE-^+EfA 

50 

Since the indenter is rigid, substitution of E*^ in Eq. (8) into 
the above equation gives 

55 

ho "UJ ~\E„r[i-a + K + 4Ka^)9 + XK€f^]) 

60 

For different E^^ and coating thickness, Eq. (10) and dis- 
placement curves in FIG. 2 are used to determine a so that 
Eq. (10) with these a values gives almost identical displace- 
ment to FIG. 2. FIG. 3 depicts the results of oJT^ against T^ 
65 with different coating/base modulus ratios E^„^. It is noted 
that in the range of small T^, value of a/T^ is strongly 
dependent on the ratio of E^^. Curves with E^^ greater than 
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one are above those with E„„. less than one. Curves with two 
extreme E^,,^. vakies (0.25 and 4) are close to unity. Some 
negative values are obtained for a when E^,.<1. It is of 
interest to find representative a which could be a function of 
coating thickness, but should be independent of the ratio 5 
E^^. If this kind of a rather than those exact values depicted 
in FIG. 3 is used to predict displacement with Eq. (10), 
discrepancy between the predicted displacement and that in 
FIG. 2 is expected to occur and should be acceptably small. 
FIG. 4 presents such tw o sets of a which are picked out from 10 
FIG. 3 w ith E^,,^.=3. Error used to quantify the discrepancy is 
defined as the difference in percentage between the values 
predicted with Eq. (10) and the values in FIG. 2 (numerical 
solutions) normalized by the latter. The relative error shown 
in FIG. 5 demonstrates the discrepancy between predicted 15 
and numerical ones is within one percent for E^^l and 
within two percent for E^^<1. In order to facilitate data 
processing, curves in FIG. 4 are replaced by the following 
simple expressions, 

20 



Tc 



1.065, Tc<025 

0.948 + 0.6377; - 0.683^,^ 0.25 < T, < 0.5 

1.064 + 0.2117; - 0.287r^^ 0.5 < < 0.9 
1 .328 - 0.382r, + 0.0466^2 0.9 <, Tc 



(11) 
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5. If necessary, a calibration sample with known material 
properties, e.g. fused silica, can be used to calibrate the 
modified Young' s moduli or the radius of the indenter. 
After the calibration sample is indented by the indenter, 
the load-displacement curve is obtained. In the calibra- 
tion with known indenter' s radius, Eq. (1) is used to 
find E^. ^, which is then used to determine E*^. from Eq. 
(2). In the calibration with known modified Young's 
moduli of the indenter, Eq. (2) is used to find E. ^, which 
is then used to detennine R from Eq. (1) and Eq. (4). 

6. Estimate the maximum load and make sure the load is 
applied within the elastic range. Record one or several 
data sets of the load and the displacement. 

7. is evaluated by Eqs. (5) and (7) with load values 
corresponding to measured data points. 

8. Use Eq. (1) to obtain E^. ^ with load and displacement 
values corresponding to measured data points. 

9. Use Eq. (2) to find E*,. 

10. Equation (8) is used to find Young's modulus of the 
coating E^, and Poisson's ratio of the coating v^, and to 
find similar properties of the base E^ and as well if 
they are unknown. Conventional minimization analysis 
can be used to find values of these properties by 
minimizing (i.e minimizing error) the following 
expression rewritten from Eq. (8), 



The corresponding relative error with this function of a is 
similar to what is depicted in FIG. 5. As mentioned before, 30 
the validity of Eq. (11) should be tested for changing the 
following conditions: Poisson's ratio of coating, load, 
indenter' s radius, and Young's modulus of the indenter. 
Several such tests have been carried out. It is found that the 
change of load and indenter' s radius has no effect on the 35 
relative error and the change of Poisson's ratio and Young's 
modulus of the indenter slightly affects the relative error. 
Therefore, the extended Hertzian analysis and Eq. (11) 
pursuant to an embodiment of the invention reflect the 
relationship between load and displacement to a satisfactory 40 
degree. 

The fraction on the right-hand side of Eq. (8) is corre- 
sponding to the presence of the substrate. Note that the 
Young's modulus expression of the sample includes the 

Young's moduli and the Poisson's ratios of the coating and ^5 
the substrate and a which is related to coating thickness. 
Similarly, if the indenter is covered by a coating (FIG. 6), 
E*. in Eq. (2) should be expressed in the same way as E*^. 
Therefore, coated indenters could be used to probe the 
sample's modulus. 

A method to determine Young's modulus and Poisson's 
ratio of the coating and the substrate pursuant to an embodi- 
ment of the invention is described as follows and illustrated 
in FIG. 13: 

1. The sample is prepared and the indenter instrument 
properly calibrated. 

2. The coating thickness of the sample is known or 
measured with proper methods. Choose a spherical 
indenter with a radius much greater than the coating 
thickness. The indenter could be made of diamond, 
sapphire, ruby, tungsten carbide, or other materials. 

3. Record the modified Young's modulus and radius of the 
indenter, if these data are available. 

4. If the radius of the indenter is not available, one could 65 
use instruments, such as profilometer or scanning elec- 
tron microscopy, to determine it. 



where a values are expressed in Eq. (11). 

11. Average Young's modulus and Poisson's ratio values 
obtained in step 10 and evaluate E^^. 

12. In order to improve the accuracy, particularly for 
E^^<1, the E^^ value is used to detennine a value from 
FIG. 3 by interpolation. 

13. Repeat steps 10 and 11. 

EXPERIMENTAL EXAMPLES 

Two samples were prepared to demonstrate practice of 
embodiments of the invention. Commercial silicon wafers 
were the substrates on which TiN was deposited with two 
different thickness. A ball with a radius of 20 mm (FIG. 7c) 
was used to wear away the coatings and part of the substrate. 
Wear scars shown in FIGS. l{a) and (/?) were used to 
determine the coating tliickness. Sample Al and A2 have 
2.08 and 3.15 micron thick coatings, respectively. Both 
samples have very smooth surfaces (Rq is around 10 nm) 
except some scattering particles (see FIG. 7 around the wear 
scars). Load and displacement curves were obtained by 
using a commercial nano-indentation system and a sapphire 
sphere with the following parameters: Young's modulus of 
441 GPa, Poisson's ratio of 0.234, and radius of 122 microns 
(E'^'. =4 6 6. 5 5 GPa). Nine tests are carried out and average 
indentation data were shown in FIG. 8 with error bar 
standing for the standard deviation. For the silicon wafers, 
Poisson's ratio is assumed as 0.27 (Reference: Chudoba et 
al. 2000) and Young's modulus is measured to be around 159 
GPa (E*^=171.50 GPa). By using these properties, coating 
thickness, and Eqs. (1) and (2), FIG. 8 is transformed into 
FIG. 9 in a form of modified Young's modulus of samples 
vs normalized coating thickness. It is noted that two samples 
have very similar behavior as expected when T^. is less than 
1, corresponding to displacement for Al greater than 23 nm. 
Although Al has significant larger E*^ than A2 at T^l, the 
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discrepancy is probably due to several uncertainties, such as 
roughness and shape of indenter tip. Material properties (E 
v^,, E^, v^) of sample Al and A2 are obtained by minimizing 
Eq. (12) with all indentation data from 23 to 200 nm 
displacement and are shown in FIG. 10a and 10b and 11a 5 
and 11^, respectively. The average Young's modulus for 
each TIN coating is 383 or 376 GPa, which agree with each 
other well and are close to the measurements obtained by 
Reference: Schneider and Tucker (1996). The average Pois- 
son's ratio for each TiN coating is 0.29. Also, the average lo 
Young's moduli of two silicon substrates measured for two 
samples are very close: 159 and 163 GPa, respectively. 

The invention is not limited to practice using the particu- 
lar substrates and coatings and indenters described in the 
above examples and can be practiced with a wide variety of 15 
substrates having a wide variety of coatings thereon using 
diflFerent indenters to determine mechanical properties 
thereof as described. 



20 



Indentation Data Probed by Other Types of Indenters 

Pursuant to another embodiment of the invention, Equa- 
tion 8 [Eq. (8)] is applicable in the elastic regime for other 
types of indenters, beside the spherical shape. The param- 
eters a in Eq. (8) is expected to be a function of non- 
dimensional coating thickness and depends on the indenter 
shape. Note that the coating thickness is normalized by 
characteristic length determined from projected area (A) of 
contact between the indenter and the sample: T^=t^yjt/A. 
With a given indenter shape, a values could be determined 
beforehand with the same procedure described above. In 
order to be simple, the a values for spherical indenters could 
be adopted here approximately. For sharp indenters, the 
uncertainty in the elastic loading portion may be high due to 
the limit of the instruments. Data in the elastic unloading 
portion could be used in our approaches with appropriate 
consideration of the impression due to plastic deformation. 
Similar to the procedure for the spherical indenter, the 
procedure used to determine Young's modulus and Pois- 
son's ratio of the coating and the substrates from the 
indentation data probed by any indenters is described as 
follows: 

1. The sample is prepared and the indenter instrument 
properly calibrated. 

2. The coating thickness of the sample is known or 
measured with proper methods. Choose an indenter 
made of diamond, sapphire, ruby, tungsten carbide, or 
other materials. 

3. Record the modified Young's modulus of the indenter, 
if it is available. 

4. A calibration sample with known material properties, 50 
e.g. fused silica, can be used to calibrate the area 
function of the indenter. After the calibration sample is 
indented by the indenter, the load-displacement curve is 
obtained. Conventional methods are used to find the 
area function. 

5. Record one or several curves of the load and the 
displacement. 

6. Use conventional methods to obtain A and E. ^ with load 
and displacement values corresponding to measured 
data points in the unloading portion. go 

7. T^ is evaluated by T^=t^^Jt/A. 

8. Use Eq. (2) to find E*,. 

9. Equation (8) is used to find E^ and v^, and to find E^ and 
v^ as well if they are unknown. Conventional minimi- 
zation analysis can be used to find values of these 65 
properties by minimizing (i.e. minimizing error) the 
following expression rewritten from Eq. (8), 



55 



12 



(12) 



Although the invention has been described in connection 
with certain embodiments thereof, those skilled in the art 
will appreciate that the invention is not so limited and that 
changes, modifications, and the like can be made thereto 
within the scope of the invention as set forth in the appended 
claims. 
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What is claimed is: 

1. A method of determining Young's modulus of a coating 
on a substrate comprising a coated sample, comprising the 
steps of generating load versus displacement indentation 25 
data using an indenter to indent the coated sample in an 

elastic range of the coating; determining a modified Young's 
modulus of the coated sample (E*J using the indentation 
data; and determining the Young's modulus of the coating 
(E^) using a relationship that expresses the modified Young' s 30 
modulus of the coated sample (E*^) as a ftinction of Young's 
modulus and Pois son's ratio of the coating. Young's modu- 
lus and Poisson's ratio of the substrate, and a parameter 
related to thickness of the coating. 

2. The method of claim 1 wherein the relationship is 35 
expressed as: 



7. Tlie method of claim 6 wherein the partial spherical 
surface has a radius greater than the thickness of the coating. 

8. The method of claim 1 wherein the indenter includes an 
uncoated surface. 

9. The method of claim 1 wherein the indenter includes a 
coated surface. 

10. In a method of determining Young's modulus of a 
coating on a substrate, the improvement comprising using 
the following relationship in determination of the Young's 
modulus: 



\-{X-\-k-^Akq^)9 + Xk9^ 
l+4ar/^0-A/t02 



(8) 



is the shear modulus ratio; 



■1 



r + (3-4v,) 



and A = 1 - 



4(1 -vc) 
l+r(3-4vc)' 



1 -(A + /f + W)6/ + A/f^ 
1 +AaKe-XK^ 



(8) 



where 



is the shear modulus ratio; 



■ 1 



r + (3-4v,) 



and A = 1 - 



4(1 -Vc) 
l+r(3-4vc)' 
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and 6=exp(-2a) where parameter a is related to non- 
dimensional coating thickness. 

3. The method of claim 1 further including determining 

the Poisson's ratio of the coating using the relationship. 

4. The method of claim 1 fiirther including determining 
Young's modulus of the substrate using the relationship. 

5. The method of claim 1 further including determining 
Poisson's ratio of the substrate using the relationship. 

6. The method of claim 1 wherein the indentation data is 
generated using the indenter having at least a partial spheri- 
cal surface. 
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and 6=exp(-2a) where parameter a is related to non- 
dimensional coating thickness, said relationship expressing 
the measured modified Young's modulus of the coated 
substrate (E*J as a function of Young's modulus and 
Poisson's ratio of the coating. Young's modulus and Pois- 
son's ratio of the substrate, and parameter a. 

11. Tlie method of claim 10 wherein parameter a is related 
to non-dimensional coating thickness T^ where T^=t^aQ and 
t^ is measured coating thickness and aO is contact radius of 
the indentation, said parameter a being determined by using 
indentation displacement data and Equation (10). 

12. In a method of determining Young's modulus of a 
coating on a substrate comprising a coated sample, the 
improvement comprising the steps of generating load versus 
displacement indentation data using an indenter to indent the 
coated sample in an elastic range of the coating; determining 
a modified Young's modulus of the coated sample (E*J that 
is dependent on coating thickness using the indentation data; 
and determining the Young's modulus of the coating itself 
using the modified Young's modulus of the coated sample 

13. The method of claim 12 wherein the indenter has at 
least a partial spherical surface. 

14. The method of claim 13 wherein the partial spherical 
surface has a radius greater than the thickness of the coating. 

15. The method of claim 12 wherein the indenter includes 
an uncoated surface. 

16. The method of claim 12 wherein the indenter includes 
a coated surface. 

17. The method of claim 12 wherein the modified Young's 
modulus of the coated sample (E* J is a function of dimen- 
sionless coating thickness. 



